Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  756 1598 (N20E34 L=134)  1131    259        4 ( 710/1543- 824/1665) -14.0 E39-E27 N15-N22
      2 1818 1202 (N29W48 L=215)   362     73        0 (1790/1174-1838/1222)   0.0 W44-W50 N28-N30
      3 1858 1111 (N25W54 L=221)   493    101        1 (1827/1082-1880/1136)   0.0 W50-W56 N23-N26
      4 1754  846 (N06W51 L=219)   396     90        0 (1718/ 821-1780/ 863)  18.4 W47-W53 N05-N08
      5 1082  796 (S17W10 L=177)   704    172        1 ( 991/ 751-1177/ 853)  51.3 W05-W13 S21-S11
      6 1765  770 (N03W55 L=223)   401     73        1 (1738/ 741-1783/ 801) -45.0 W53-W57 N01-N04
      7 1729  614 (S06W60 L=227)  1182    141        0 (1680/ 553-1760/ 685)  90.0 W56-W64 S10-S01
      8 1271  601 (S22W28 L=195)   680    256        1 (1178/ 553-1380/ 675)   0.0 W19-W36 S22-S19
      9 1454  630 (S14W38 L=205)   266     74        0 (1441/ 597-1476/ 660) -63.4 W35-W39 S16-S12
     10  807  495 (S43W03 L=171)   875    252        0 ( 715/ 469- 945/ 540)  19.4 E04-W13 S45-S38
     11 1115  413 (S38W28 L=195)   247     51        1 (1086/ 400-1149/ 431)  45.0 W26-W29 S39-S36
     12 1481  314 (S30W58 L=225)  1426    302        4 (1402/ 287-1574/ 354)  17.4 W49-W65 S32-S25
  

    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