Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1768 (N39E35 L=297)   895    153        2 ( 878/1700- 935/1831) -34.7 E42-E28 N34-N43
      2  630 1678 (N25E48 L=284)   685    139        3 ( 595/1624- 663/1743) -26.6 E56-E42 N21-N28
      3 1462 1716 (N56W12 L=344)  1179    112        2 (1408/1688-1521/1738)  37.6 W06-W20 N50-N60
      4  158 1008 (S25E67 L=264)  1380    182        4 ( 115/ 922- 216/1072)  39.8 E77-E58 S30-S20
      5  489  768 (S33E29 L=303)   694    192        2 ( 408/ 741- 589/ 790)  21.8 E36-E21 S35-S28
      6  998  623 (S28W10 L=341)   598    138        1 ( 970/ 554-1033/ 684) -83.7 W08-W11 S32-S23
      7 1324  489 (S26W36 L=  8)   823    325        1 (1202/ 442-1401/ 583)  59.7 W28-W40 S31-S18
      8  401  438 (S61E45 L=286)  1182    115        3 ( 351/ 405- 448/ 490)  -3.2 E55-E36 S63-S55
      9  503  386 (S62E29 L=303)   775     69        2 ( 477/ 348- 524/ 420) -81.9 E33-E24 S65-S58
     10  686  353 (S57E03 L=328)   802    152        1 ( 624/ 324- 745/ 403)  39.3 E10-W01 S61-S51
     11  867  236 (S60W23 L=355)   399     69        0 ( 849/ 214- 893/ 268) -21.8 W18-W28 S62-S57
  

    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