Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1211 (N12W24 L=  3)   604    200        2 (1336/1140-1427/1292) -78.7 W22-W26 N07-N17
      2  251 1165 (N22E57 L=281)   351     37        1 ( 230/1142- 275/1186) -26.6 E59-E54 N21-N23
      3 1621  889 (S11W37 L= 15)   343     72        1 (1593/ 849-1644/ 925)  53.1 W34-W39 S13-S09
      4 1853  754 (S26W63 L= 41)   645     72        0 (1825/ 735-1876/ 792)  31.0 W58-W66 S27-S23
      5 1527  656 (S25W29 L=  7)   706    213        1 (1471/ 576-1584/ 746)  48.0 W24-W33 S29-S19
      6 1358  644 (S22W16 L=355)   247     80        0 (1337/ 608-1375/ 684)  59.0 W14-W17 S24-S19
      7 1034  466 (S30E08 L=330)   383     96        0 (1016/ 414-1064/ 506)  56.3 E10-E05 S33-S26
      8  377  425 (S26E65 L=273)  1335    125        3 ( 362/ 367- 411/ 481)  25.5 E79-E58 S31-S21
      9  794  293 (S40E32 L=306)   379     58        0 ( 767/ 271- 826/ 314)   9.5 E35-E29 S41-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