Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1543 1147 (N16W31 L=284)  4769    928        8 (1244/ 780-1754/1454) -32.0 W06-W55 S02-N29
      2  741 1064 (S04E18 L=236)   227     67        0 ( 721/1031- 762/1100) -59.0 E18-E15 S05-S01
      3  337 1010 (S14E46 L=208)   420     52        0 ( 314/ 981- 358/1043)  76.0 E47-E43 S15-S11
      4  359  886 (S21E42 L=211)   425    131        0 ( 314/ 847- 401/ 933)  45.0 E46-E39 S24-S17
      5  496  860 (S20E31 L=222)   234     84        0 ( 483/ 832- 515/ 910) -63.4 E32-E29 S22-S17
      6 1610  732 (S08W44 L=298)   782    194        1 (1537/ 708-1691/ 762)  18.4 W38-W50 S10-S05
      7  763  702 (S25E11 L=243)   458    123        1 ( 699/ 675- 823/ 721)   0.0 E15-E06 S26-S23
      8 1295  628 (S20W25 L=279)   565    181        0 (1213/ 608-1377/ 665)   4.8 W18-W30 S20-S17
      9 1306  279 (S43W43 L=297)   425     53        1 (1279/ 260-1335/ 298)  33.7 W41-W44 S44-S41
  

    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