Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1371 1689 (N45W01 L= 13)  2349    341        2 (1261/1563-1493/1796) -26.6 E10-W14 N39-N52
      2  444 1707 (N19E67 L=305)   751     87        2 ( 402/1672- 472/1738) -80.5 E70-E63 N16-N22
      3  776 1652 (N23E35 L=337)  1228    288        6 ( 690/1542- 876/1726)  -9.5 E42-E24 N17-N27
      4 1179 1599 (N31E08 L=  4)   495    121        0 (1153/1540-1195/1654) -54.5 E11-E05 N27-N34
      5 1435 1444 (N29W13 L= 25)   275    118        0 (1385/1402-1473/1474) -14.0 W08-W16 N27-N30
      6 1864 1178 (N29W56 L= 67)   551     90        1 (1839/1143-1892/1220) -26.6 W51-W60 N27-N30
      7  178 1099 (S22E64 L=307)  1431    154        6 ( 124/1039- 210/1159)  32.5 E73-E59 S25-S18
      8  296  996 (S26E49 L=323)   361    103        0 ( 255/ 961- 334/1034)  40.6 E52-E45 S29-S23
      9  807  861 (S22E08 L=  4)   288     76        1 ( 760/ 843- 852/ 873)  21.8 E10-E05 S23-S20
     10  227  790 (S40E58 L=314)   616    116        0 ( 186/ 748- 276/ 818)  17.4 E66-E50 S43-S38
     11 1712  647 (S06W56 L= 68)   324     36        1 (1694/ 625-1728/ 669) -26.6 W54-W58 S07-S04
     12 1427  402 (S28W50 L= 61)  1357    261        0 (1358/ 309-1534/ 470)  85.6 W46-W54 S35-S21
  

    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