Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1605 (N37E13 L=355)   368    144        2 ( 931/1560-1032/1651) -26.6 E17-E07 N34-N40
      2  546 1362 (N15E36 L=331)   726    297        1 ( 505/1225- 571/1477) -53.7 E43-E32 N06-N21
      3  262 1347 (N10E62 L=305)  1182    169        4 ( 238/1268- 295/1424) -68.2 E66-E58 N04-N14
      4 1701 1331 (N28W47 L= 54)   558    105        0 (1658/1287-1736/1372) -33.7 W42-W50 N25-N30
      5  798 1198 (N08E16 L=351)  1004    256        1 ( 746/1084- 847/1303)  81.9 E17-E14 N00-N14
      6 1025 1169 (N09E02 L=  6)   269    103        0 ( 999/1123-1052/1215) -59.0 E03-W00 N06-N11
      7  481 1104 (S02E36 L=332)   334     70        0 ( 451/1079- 517/1134) -21.8 E38-E33 S02-N00
      8 1020  654 (S22W05 L= 12)  2202    673        1 ( 761/ 555-1295/ 772)  32.8 E11-W20 S32-S11
      9  339  625 (S35E50 L=318)   538    118        0 ( 316/ 579- 372/ 679) -66.8 E51-E46 S37-S30
     10  526  620 (S32E31 L=337)   427     93        0 ( 478/ 604- 572/ 634)   0.0 E34-E26 S33-S30
     11  994  183 (S62W19 L= 27)   425     26        0 ( 972/ 162-1013/ 203) -24.0 W15-W24 S64-S59
  

    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