Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1260 (N16W05 L=284)  3430    967        7 ( 874/ 881-1452/1591) -33.1 E21-W28 S01-N31
      2 1476 1547 (N42W23 L=302)   327     79        0 (1448/1523-1509/1561)  21.8 W20-W25 N39-N43
      3  288 1344 (N04E59 L=221)   391     97        1 ( 265/1320- 316/1368) -21.8 E61-E56 N03-N05
      4  183 1079 (S13E63 L=217)   546    118        2 ( 165/1030- 204/1134)  81.9 E64-E60 S16-S09
      5  389 1004 (S14E41 L=238)   328     66        2 ( 367/ 959- 408/1040)  90.0 E42-E39 S16-S11
      6  436  817 (S25E35 L=244)   644    204        1 ( 357/ 786- 507/ 849)   0.0 E42-E29 S26-S22
      7  971  779 (S17W01 L=281)  1699    536        5 ( 733/ 732-1225/ 844)  24.9 E13-W15 S21-S08
      8 1411  356 (S35W46 L=325)  4880    738       13 (1129/ 312-1642/ 437)  26.6 W25-W62 S43-S24
      9 1056  357 (S43W20 L=299)   325     38        0 (1026/ 341-1079/ 377)  23.4 W17-W22 S45-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