Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1335 1655 (N45W02 L=169)   380     99        0 (1284/1639-1380/1668)  23.2 E01-W06 N42-N46
      2  461 1548 (N13E53 L=113)   352     77        2 ( 434/1511- 488/1572) -20.6 E56-E48 N11-N14
      3 1407 1550 (N39W11 L=177)   316     69        0 (1372/1530-1439/1570)  -9.5 W07-W13 N37-N40
      4 1646 1303 (N30W34 L=200)  1962    370        3 (1500/1199-1782/1414)  -8.1 W20-W48 N26-N32
      5 1060 1087 (N01W00 L=167)   415    136        0 (1022/1036-1089/1147) -33.7 E03-W03 S01-N03
      6  173 1082 (S20E65 L=101)   728    130        0 ( 146/1020- 209/1133)  48.8 E69-E61 S24-S16
      7 1628 1103 (N17W35 L=202)   285     64        2 (1608/1076-1649/1128) -26.6 W33-W37 N15-N18
      8 1058  903 (S10W05 L=171)   851    279        1 ( 970/ 809-1131/1001)  75.1 W02-W07 S17-S02
      9  319  890 (S29E46 L=121)   328     73        1 ( 304/ 854- 340/ 921) -63.4 E47-E43 S30-S26
     10  176  829 (S36E66 L=100)   599     98        0 ( 163/ 782- 187/ 881)  45.0 E70-E62 S38-S32
     11 1448  560 (S18W40 L=207)   259     25        0 (1431/ 541-1469/ 579)  90.0 W38-W41 S19-S17
  

    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