Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  762 1725 (N33E41 L=268)  1013    139        0 ( 709/1681- 816/1781) -14.0 E50-E34 N30-N36
      2  549 1630 (N22E51 L=258)   798    156        1 ( 512/1564- 572/1695) -31.6 E59-E46 N17-N25
      3 1368 1646 (N46W07 L=316)   635    121        0 (1328/1612-1423/1669)   0.0 W02-W12 N43-N47
      4  600 1555 (N18E42 L=267)   493    110        0 ( 562/1506- 633/1601)  90.0 E43-E40 N14-N21
      5  421 1558 (N14E57 L=252)   442     37        2 ( 396/1536- 448/1581)  21.8 E59-E54 N13-N16
      6 1744 1215 (N27W45 L=354)   794     87        1 (1713/1171-1764/1264) -59.0 W41-W47 N24-N29
      7 1280 1193 (N13W11 L=320)   314     95        0 (1258/1156-1307/1248) -45.0 W08-W12 N10-N15
      8 1843 1158 (N26W57 L=  6)   890    145        1 (1801/1098-1883/1219) -22.6 W51-W63 N23-N29
      9  531 1090 (S11E32 L=277)   420    103        0 ( 488/1069- 582/1109)   8.1 E34-E27 S12-S09
     10  413  762 (S33E36 L=273)   406     75        0 ( 382/ 722- 438/ 801)  56.3 E39-E34 S36-S30
     11 1256  590 (S22W26 L=335)   363    122        1 (1237/ 535-1284/ 646) -35.5 W23-W30 S24-S19
     12  612  555 (S42E17 L=292)  1260    260        4 ( 523/ 507- 741/ 591)  20.2 E25-E06 S47-S37
  

    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