Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  527 1618 (N23E51 L= 42)   533    105        2 ( 478/1594- 582/1639)   0.0 E56-E45 N21-N24
      2 1546 1490 (N39W27 L=120)   388    116        1 (1505/1457-1600/1510)  -6.3 W23-W32 N37-N40
      3 1011 1341 (N16E08 L= 85)   583    219        2 ( 956/1279-1116/1384)  11.3 E11-E01 N11-N19
      4 1293 1318 (N21W10 L=103)   243     53        1 (1273/1294-1310/1346)  76.0 W08-W11 N19-N23
      5 1290 1245 (N17W11 L=104)   498    141        2 (1244/1202-1335/1270)  31.0 W08-W14 N13-N18
      6  941  609 (S29W04 L= 97)  1209    437        2 ( 761/ 543-1113/ 710)  45.0 E06-W12 S37-S18
      7  804  364 (S50W02 L= 95)   382     27        0 ( 779/ 342- 832/ 385) -14.0 W00-W04 S51-S48
      8 1089  302 (S46W28 L=121)  1546    301        6 ( 982/ 278-1203/ 348)  23.6 W19-W37 S50-S41
  

    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