Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1344 1625 (N42W31 L=331)   285     62        1 (1322/1599-1366/1645)  31.0 W29-W34 N39-N43
      2 1351 1450 (N29W26 L=326)   528    156        2 (1274/1424-1424/1468) -10.3 W20-W31 N27-N30
      3 1051 1252 (N19W03 L=303)   237    102        0 (1023/1189-1090/1296) -60.3 W01-W05 N14-N21
      4 1673  765 (S16W44 L=344)  1564    421        0 (1512/ 611-1785/ 950)  42.3 W31-W53 S25-S05
      5  274  845 (S04E55 L=245)   952    154        1 ( 244/ 785- 310/ 934) -63.4 E57-E51 S08-N01
      6  816  912 (S01E13 L=286)   210     84        0 ( 777/ 898- 864/ 925)  11.3 E15-E10 S01-N00
      7 1796  704 (S22W57 L=357)   594    108        1 (1763/ 667-1825/ 754)  63.4 W53-W59 S24-S18
      8 1692  553 (S32W49 L=349)   424     83        1 (1671/ 519-1708/ 585)  78.7 W47-W51 S34-S29
      9  628  523 (S26E31 L=269)   305     70        1 ( 588/ 494- 660/ 551) -36.9 E33-E29 S28-S23
  

    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