Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  985 1905 (N63W49 L=199)   548     58        0 ( 965/1873-1006/1937)   8.1 W40-W61 N60-N65
      2  424 1459 (N47E34 L=116)   420     45        1 ( 402/1432- 450/1482)  33.7 E36-E31 N44-N48
      3 1008 1322 (N22W07 L=157)   212     78        0 ( 970/1304-1049/1337) -11.3 W04-W09 N20-N22
      4  835 1291 (N25E04 L=146)   236     75        1 ( 822/1248- 858/1329)  63.4 E05-E02 N22-N26
      5  341 1116 (N26E43 L=107)   688    214        0 ( 297/1041- 421/1158)  10.3 E47-E36 N22-N28
      6 1688  902 (S20W38 L=188)   492    163        0 (1644/ 841-1748/ 956)  21.8 W33-W44 S22-S18
      7 1622  685 (S32W30 L=180)   502    151        0 (1561/ 645-1680/ 741)  15.3 W24-W35 S32-S29
      8  833  563 (S17E24 L=127)   269     98        0 ( 789/ 535- 872/ 598)   9.5 E26-E20 S17-S15
      9  433  525 (S10E53 L= 97)   814    156        1 ( 413/ 446- 461/ 590)  56.3 E56-E48 S15-S06
     10 1409  490 (S38W10 L=160)   714    154        0 (1341/ 449-1479/ 538)   4.8 W04-W16 S39-S36
     11  574  288 (S28E59 L= 92)  1555    208        2 ( 503/ 245- 655/ 335) -12.1 E65-E47 S31-S23
  

    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