Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1446 1474 (N35W20 L=178)  3106    642        3 (1158/1344-1693/1621)  -6.3 E04-W41 N30-N39
      2  734 1553 (N25E31 L=127)   469    127        2 ( 694/1508- 772/1603) -32.0 E35-E27 N22-N27
      3 1155 1255 (N15W03 L=161)   239     67        0 (1144/1218-1169/1291) -63.4 W02-W04 N12-N16
      4 1198 1142 (N09W08 L=166)   222     94        0 (1164/1122-1253/1172)   0.0 W05-W11 N07-N10
      5  925 1065 (S01E07 L=151)   385     78        1 ( 880/1035- 961/1099)  36.9 E08-E04 S03-N00
      6  342  768 (S30E43 L=115)   409     72        0 ( 321/ 734- 364/ 799)  51.3 E45-E41 S32-S27
      7  951  772 (S18W01 L=159)   215     23        1 ( 932/ 752- 973/ 793)  71.6 W00-W01 S19-S16
      8 1490  726 (S10W35 L=193)  1344    279        4 (1362/ 692-1593/ 779)  14.0 W26-W42 S12-S05
      9 1737  674 (S08W58 L=215)   710    123        1 (1697/ 620-1763/ 718) -45.0 W54-W60 S10-S05
     10 1434  413 (S30W44 L=202)   329     83        1 (1395/ 399-1465/ 428)  11.3 W41-W46 S31-S28
     11 1330  256 (S44W48 L=206)   399     67        0 (1315/ 233-1345/ 284) -38.7 W45-W50 S45-S41
  

    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