Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1001 1367 (N11E10 L= 13)   373    117        1 ( 959/1308-1038/1413) -35.5 E13-E06 N08-N13
      2  305 1110 (S18E49 L=334)  1321    293        2 ( 248/ 981- 346/1265)  84.0 E52-E46 S26-S07
      3  644 1113 (S11E24 L=359)   807    289        7 ( 538/1062- 771/1157)   3.6 E31-E15 S13-S09
      4 1878 1011 (N17W58 L= 81)   825    118        3 (1856/ 949-1895/1062) -74.1 W55-W60 N13-N20
      5 1059  803 (S19W08 L= 31)   498    233        0 (1011/ 706-1094/ 911) -54.0 W02-W11 S24-S13
      6 1653  793 (S03W45 L= 67)  1249    267        2 (1576/ 747-1745/ 871)  -6.7 W36-W53 S04-N00
      7  796  719 (S30E06 L= 17)  1123    418        2 ( 599/ 679-1001/ 756)  17.1 E19-W07 S34-S25
      8 1245  706 (S19W21 L= 44)   211     72        1 (1204/ 681-1277/ 735)  45.0 W19-W22 S21-S16
      9 1702  655 (S08W54 L= 77)   342     56        0 (1681/ 621-1721/ 686)  90.0 W53-W55 S10-S06
     10 1256  535 (S28W29 L= 52)   536    213        1 (1175/ 517-1379/ 558)  30.3 W23-W35 S30-S23
     11 1507  454 (S25W49 L= 72)  2670    530        4 (1328/ 385-1713/ 549)  42.3 W39-W62 S33-S13
     12 1612  398 (S24W61 L= 84)  1057    208        6 (1579/ 338-1644/ 467) -53.1 W57-W64 S27-S19
  

    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