Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1428 1469 (N24W38 L=199)  1637    292        3 (1305/1374-1508/1568) -72.6 W33-W41 N17-N33
      2 1626 1510 (N22W56 L=217)   684    195        2 (1536/1473-1695/1530) -24.8 W47-W61 N17-N24
      3 1291 1183 (N11W19 L=180)   578    331        0 (1233/1121-1380/1228)   0.0 W14-W25 N06-N13
      4 1398  948 (S05W21 L=182)   395    142        2 (1329/ 929-1463/ 969) -35.5 W16-W24 S07-S02
      5 1669  749 (S24W38 L=199)  1847    480        3 (1528/ 600-1806/ 931)  28.3 W25-W53 S30-S15
      6 1521  857 (S13W28 L=189)   345    101        0 (1498/ 815-1536/ 903)  59.0 W25-W29 S15-S10
      7  455  870 (N07E39 L=122)   367     87        0 ( 416/ 852- 496/ 888)  -9.5 E42-E36 N05-N07
      8 1752  697 (S30W47 L=208)  1130    215        1 (1682/ 604-1816/ 772)  29.7 W40-W54 S34-S26
      9  649  580 (S14E33 L=128)   291     77        1 ( 610/ 563- 687/ 603)   0.0 E35-E29 S14-S12
     10  353  499 (S15E62 L= 99)   456     58        2 ( 340/ 470- 370/ 528)  29.7 E65-E58 S17-S13
     11  568  474 (S19E43 L=118)   289     67        1 ( 548/ 441- 588/ 499)  76.0 E45-E40 S21-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