Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  897 1780 (N60W11 L=287)   370     36        1 ( 876/1754- 915/1805)  39.8 W08-W15 N57-N62
      2 1439 1605 (N33W45 L=321)  1925    277        4 (1336/1563-1573/1669) -42.7 W39-W54 N27-N39
      3  595 1459 (N41E24 L=252)   996    241        0 ( 480/1441- 714/1482) -12.8 E35-E13 N37-N43
      4  241 1456 (N45E71 L=205)  1069     90        1 ( 205/1428- 292/1475)   2.6 E84-E60 N43-N46
      5 1743 1326 (N09W56 L=332)   357     85        0 (1718/1281-1763/1365)  23.2 W52-W59 N07-N10
      6 1737  668 (S30W46 L=321)  2793    395        4 (1610/ 565-1850/ 831)  20.6 W33-W57 S36-S22
      7  216  634 (S07E71 L=205)  1078     95        5 ( 191/ 582- 238/ 673)  56.3 E75-E66 S10-S04
      8 1459  642 (S25W21 L=297)   432     89        0 (1416/ 624-1502/ 669)   0.0 W17-W24 S26-S23
      9  266  580 (S10E67 L=209)   585     78        2 ( 245/ 544- 283/ 615)  45.0 E69-E64 S12-S07
     10  632  448 (S23E38 L=237)   437     63        0 ( 601/ 419- 669/ 480) -63.4 E40-E36 S24-S20
     11 1039  269 (S44E17 L=259)  3380    500        6 ( 846/ 186-1260/ 389)   9.5 E39-W03 S48-S38
  

    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