Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  495 1605 (N40E49 L=353)   448     76        1 ( 464/1579- 526/1630)  59.0 E51-E46 N37-N42
      2  961 1454 (N31E05 L= 38)   296     69        2 ( 928/1433- 998/1468)  21.8 E07-E02 N29-N31
      3  521 1225 (N15E34 L=  9)   836    242        3 ( 481/1117- 580/1321) -53.1 E38-E28 N08-N20
      4  182 1262 (N15E71 L=332)  1024    111        1 ( 156/1225- 211/1317) -31.0 E78-E64 N12-N18
      5 1274 1254 (N18W16 L= 59)   478    132        1 (1224/1223-1334/1276)  14.0 W12-W20 N15-N19
      6 1829  837 (S09W62 L=104)   804    137        0 (1794/ 777-1871/ 892)  45.0 W58-W66 S12-S05
      7  532  835 (S10E32 L= 11)   391     76        1 ( 495/ 806- 567/ 873)  31.0 E34-E29 S11-S07
      8  723  566 (S27E21 L= 22)   308    100        1 ( 707/ 520- 739/ 618)  81.9 E21-E19 S30-S23
      9 1527  437 (S36W43 L= 86)  1311    192        3 (1488/ 346-1561/ 521) -66.8 W38-W49 S43-S29
  

    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