Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1836 (N57E51 L=309)  1421     84        2 ( 684/1802- 771/1865)  33.7 E59-E45 N53-N60
      2  853 1617 (N40E19 L=340)   937    228        3 ( 760/1568- 959/1658) -11.9 E28-E09 N36-N42
      3  396 1524 (N28E56 L=303)   669    120        1 ( 353/1488- 442/1568) -14.9 E64-E49 N26-N31
      4 1334 1436 (N30W19 L= 19)   260     93        1 (1311/1402-1363/1478) -59.0 W17-W21 N28-N33
      5  237 1418 (N20E71 L=289)  1564    103        4 ( 207/1371- 269/1464) -16.7 E86-E63 N16-N22
      6 1112 1443 (N29W03 L=  3)   237     53        1 (1086/1428-1139/1463) -33.7 W01-W04 N28-N30
      7  634 1331 (N18E28 L=331)   601    283        0 ( 537/1240- 730/1400) -25.2 E37-E20 N13-N21
      8 1717 1178 (N16W48 L= 48)   331     58        1 (1702/1138-1731/1207) -63.4 W47-W50 N13-N17
      9 1153 1080 (N06W08 L=  7)   346     79        1 (1130/1044-1181/1109)  53.1 W06-W09 N04-N08
     10 1710  645 (S18W55 L= 54)   475     48        1 (1687/ 628-1738/ 668)   0.0 W52-W58 S18-S16
     11 1473  556 (S25W36 L= 36)   364    101        0 (1428/ 541-1522/ 581) -12.5 W31-W40 S25-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