Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1634 (N50W03 L= 63)   267     62        2 ( 806/1613- 863/1657)  11.3 W00-W05 N48-N51
      2 1271 1335 (N20W22 L= 82)   251     73        1 (1237/1312-1302/1354) -45.0 W20-W24 N17-N21
      3 1810 1150 (S05W57 L=117)   902    126        2 (1784/1095-1829/1197)  51.3 W53-W59 S07-S01
      4  683 1098 (N18E19 L= 41)   313     67        3 ( 654/1067- 715/1124)  14.0 E21-E17 N16-N19
      5  246  903 (N13E58 L=  2)  1799    260        4 ( 204/ 813- 316/1010) -72.9 E62-E51 N06-N19
      6  471  609 (S08E45 L= 15)   329     88        0 ( 438/ 575- 503/ 628) -45.0 E47-E43 S10-S06
      7 1214  438 (S34E01 L= 59)  1687    328        6 (1123/ 324-1304/ 583)  43.3 E10-W07 S41-S25
      8  900  432 (S27E22 L= 38)   622    175        0 ( 848/ 382- 930/ 505)  37.9 E26-E17 S29-S22
      9  533  397 (S23E51 L=  9)  1045    194        2 ( 501/ 308- 552/ 486)  47.3 E58-E45 S29-S16
     10 1032  232 (S47E23 L= 37)  1543    159        4 ( 976/ 173-1086/ 298)  30.5 E33-E16 S52-S40
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y