Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  362 1551 (N16E57 L= 33)   995    208        4 ( 309/1495- 392/1594) -60.3 E62-E52 N12-N19
      2 1712 1528 (N40W43 L=133)   423    100        0 (1690/1499-1740/1549)  45.0 W40-W46 N37-N43
      3 1528 1427 (N27W26 L=115)   873    199        2 (1430/1368-1608/1486) -17.1 W18-W31 N24-N29
      4 1469 1180 (N11W24 L=114)   572    108        1 (1439/1110-1488/1229) -66.8 W21-W26 N06-N13
      5 1689 1073 (N09W41 L=131)   608    179        3 (1637/ 994-1737/1122) -56.3 W36-W44 N04-N11
      6  264 1058 (S15E52 L= 38)   389     88        1 ( 232/1034- 300/1084)  -8.1 E55-E48 S16-S13
      7 1855  891 (N03W60 L=150)   710    182        0 (1815/ 845-1908/ 925)  29.1 W57-W66 N00-N06
      8 1066  750 (S22W08 L= 97)   468    145        0 (1003/ 717-1140/ 817)  53.1 W04-W10 S24-S16
      9 1724  566 (S18W59 L=149)  1064    161        2 (1666/ 496-1785/ 631)  60.9 W55-W63 S22-S13
     10 1066  369 (S46W19 L=109)   605     99        2 (1019/ 343-1124/ 397)   0.0 W15-W25 S47-S44
     11  862  208 (S66W11 L=101)   493     93        2 ( 839/ 180- 884/ 236) -29.1 W06-W16 S69-S63
  

    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