Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  550 1610 (N23E49 L= 41)   432    117        0 ( 505/1587- 608/1630)   0.0 E53-E42 N22-N24
      2 1621 1455 (N39W35 L=125)   456     95        0 (1580/1421-1661/1481) -14.0 W30-W38 N37-N40
      3 1352 1392 (N27W13 L=103)   261    148        1 (1289/1358-1404/1412)   7.1 W08-W16 N25-N28
      4 1086 1334 (N18E03 L= 87)   588    207        3 ( 999/1292-1162/1377)  18.4 E08-W01 N15-N19
      5 1333 1231 (N17W14 L=104)   500    132        1 (1287/1203-1379/1255)   9.5 W11-W17 N14-N18
      6  978  593 (S29W07 L= 97)  1336    441        3 ( 790/ 526-1160/ 690)  43.5 E04-W16 S36-S18
      7  834  356 (S50W04 L= 94)   332     44        0 ( 811/ 337- 861/ 377) -11.3 W02-W07 S51-S48
      8 1113  287 (S47W31 L=121)  1703    219        1 (1003/ 262-1224/ 334)  14.7 W22-W41 S50-S43
  

    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