Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1526 1471 (N20W45 L=338)   298     63        0 (1511/1438-1541/1501)  76.0 W43-W46 N18-N22
      2  935 1467 (N35W05 L=298)  1147    308        1 ( 807/1439-1084/1496) -14.7 E04-W15 N31-N38
      3 1288 1425 (N24W26 L=320)   812    299        0 (1169/1329-1379/1468) -56.3 W19-W29 N16-N28
      4 1267 1264 (N15W20 L=313)   607    198        0 (1227/1199-1322/1305)  45.0 W17-W23 N11-N18
      5 1808 1229 (S01W59 L=352)   430     60        2 (1780/1205-1834/1254)   0.0 W55-W63 S02-N00
      6  202 1135 (N28E62 L=231)   472     98        0 ( 190/1087- 218/1181)  80.5 E63-E60 N24-N30
      7 1785 1103 (S07W53 L=346)   503    106        2 (1756/1070-1808/1155)  53.1 W49-W55 S09-S04
      8  640  944 (N10E25 L=268)  1035    356        1 ( 604/ 796- 687/1104)  87.0 E26-E21 N00-N19
      9  217 1013 (N20E60 L=234)   804    133        0 ( 187/ 963- 267/1059) -37.9 E63-E54 N16-N23
     10  144  969 (N18E71 L=222)   766     69        3 ( 124/ 934- 169/1001) -63.4 E74-E66 N15-N19
     11  347  784 (N05E50 L=243)   328     96        1 ( 321/ 750- 387/ 820) -51.3 E52-E47 N02-N07
     12  278  599 (S06E65 L=228)  6050    452       12 ( 197/ 434- 371/ 776)  74.7 E81-E50 S17-N05
     13 1671  748 (S25W37 L=330)   368     69        1 (1647/ 723-1705/ 772)   0.0 W34-W40 S26-S23
  

    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