Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1036 1783 (N41E28 L=221)  1530    247        2 ( 982/1669-1117/1858) -70.7 E32-E21 N31-N51
      2 1375 1499 (N32W07 L=256)   948    165        1 (1348/1414-1398/1576) -59.0 W04-W10 N27-N37
      3 1612 1386 (N34W27 L=276)  1530    348        0 (1529/1243-1697/1513) -31.4 W18-W36 N26-N40
      4 1165 1403 (N19E03 L=246)   245     85        0 (1147/1360-1195/1434)  90.0 E03-E01 N16-N21
      5 1730 1097 (N20W42 L=291)   369     71        1 (1710/1058-1744/1130) -53.1 W40-W44 N17-N21
      6  956  812 (S19W02 L=251)   196     69        0 ( 931/ 783- 982/ 841) -14.0 W00-W04 S20-S18
      7 1038  794 (S18W07 L=256)   217     83        0 ( 994/ 781-1072/ 819)  53.1 W05-W08 S19-S15
      8  985  445 (S40W18 L=267)   547    128        1 ( 943/ 404-1014/ 491) -82.9 W15-W20 S44-S36
      9  840  406 (S48W11 L=259)   306     60        0 ( 814/ 388- 874/ 434)  -9.5 W07-W13 S49-S46
     10 1216  300 (S41W44 L=293)   527     80        0 (1172/ 273-1258/ 325)   6.3 W40-W49 S42-S39
     11 1476  289 (S31W66 L=315)   848     69        2 (1437/ 265-1513/ 311)   4.1 W60-W74 S32-S29
  

    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