Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  689 1686 (N58E17 L=  1)  1744    228        4 ( 614/1617- 798/1727)  -2.5 E26-E03 N52-N61
      2 1232 1649 (N41W32 L= 50)  1985    412        6 (1039/1617-1436/1675) -21.2 W16-W48 N35-N47
      3  271 1344 (N37E59 L=320)   556     93        0 ( 241/1309- 308/1369)  33.7 E62-E53 N34-N38
      4  490 1315 (N33E34 L=345)   313     64        1 ( 453/1294- 525/1331)  -9.5 E36-E30 N31-N33
      5 1187 1237 (N16W14 L= 32)   862    202        3 (1113/1161-1238/1305) -63.4 W10-W15 N11-N21
      6  627 1181 (N22E23 L=355)   403    112        1 ( 613/1123- 648/1241)  81.9 E24-E21 N18-N25
      7 1607  695 (S25W34 L= 52)  1305    192        5 (1577/ 616-1658/ 785) -65.2 W31-W39 S31-S18
      8 1303  690 (S19W12 L= 30)   582    147        0 (1258/ 625-1359/ 746)  40.6 W07-W15 S22-S15
  

    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