Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  793 1701 (N33E29 L=254)  1062    238        1 ( 728/1598- 823/1783) -40.9 E39-E23 N26-N39
      2  482 1634 (N25E51 L=233)   603     74        2 ( 437/1617- 526/1660)   6.3 E55-E46 N24-N26
      3 1339 1634 (N38W12 L=295)   258     74        0 (1304/1609-1374/1658) -18.4 W09-W15 N36-N39
      4 1258 1596 (N33W06 L=290)   263     96        0 (1232/1547-1277/1639) -56.3 W04-W08 N30-N36
      5  382 1538 (N18E55 L=228)  1023    136        4 ( 349/1474- 409/1593) -58.0 E60-E51 N13-N21
      6  184  987 (S17E62 L=221)  3361    357        4 ( 121/ 853- 277/1167)  57.7 E72-E52 S25-S06
      7  587  838 (S23E25 L=258)  1549    521        4 ( 407/ 742- 797/ 986)  35.1 E40-E12 S30-S11
      8  888  826 (S20E05 L=278)   561    273        1 ( 827/ 759- 973/ 871)  12.5 E09-W00 S23-S17
      9 1454  767 (S14W31 L=314)   277     45        0 (1440/ 739-1471/ 803) -45.0 W29-W32 S15-S12
     10 1544  582 (S23W44 L=327)  1527    415        0 (1392/ 475-1678/ 725)  61.2 W36-W49 S32-S12
     11  749  283 (S61E11 L=272)  1120    148        1 ( 667/ 253- 814/ 310)   6.0 E22-E03 S63-S58
  

    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