Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  857 1738 (N57W14 L= 39)   463     82        1 ( 831/1715- 900/1760) -51.3 W10-W17 N54-N59
      2 1586 1677 (N27W67 L= 92)   467     48        2 (1557/1662-1616/1692) -26.6 W63-W71 N25-N28
      3 1504 1592 (N26W52 L= 77)   766    182        1 (1431/1549-1573/1626) -45.0 W46-W55 N21-N29
      4  521 1558 (N54E28 L=358)   386     55        1 ( 485/1538- 546/1582) -14.0 E32-E24 N51-N55
      5 1415 1467 (N22W37 L= 63)   254     39        0 (1392/1444-1434/1482) -63.4 W35-W39 N20-N23
      6 1325 1428 (N22W30 L= 55)   472    132        0 (1266/1402-1381/1446) -23.2 W25-W33 N19-N24
      7 1115 1248 (N18W11 L= 36)   890    230        2 (1008/1180-1190/1313) -61.4 W06-W14 N12-N24
      8  409 1137 (N27E39 L=346)   255     48        1 ( 387/1115- 437/1152)   0.0 E40-E36 N25-N28
      9 1806 1113 (S09W55 L= 80)   346     47        0 (1787/1088-1824/1136) -76.0 W52-W56 S11-S07
     10  839 1089 (N15E09 L= 16)   278     79        1 ( 817/1050- 859/1126)  45.0 E10-E07 N13-N17
     11 1870  906 (S24W59 L= 84)  3139    493        3 (1823/ 797-1930/1047)  75.1 W52-W70 S31-S16
     12  430  994 (N18E39 L=346)   313     72        1 ( 417/ 955- 443/1035) -90.0 E39-E37 N15-N20
     13 1149  932 (S01W05 L= 30)   949    219        4 (1097/ 835-1205/1021)  37.9 W00-W10 S06-N02
     14 1317  689 (S20W09 L= 34)   429    140        0 (1249/ 651-1364/ 711) -35.5 W04-W11 S23-S17
     15 1221  514 (S28E01 L= 24)   634    189        0 (1141/ 482-1308/ 558)  -5.2 E06-W05 S29-S26
     16 1510  518 (S37W20 L= 45)   729     87        1 (1474/ 474-1544/ 553)  38.7 W16-W22 S40-S34
     17 1374  399 (S42W07 L= 32)   402     62        1 (1351/ 376-1407/ 424)  21.8 W04-W09 S43-S39
     18 1276  351 (S43E02 L= 23)   345     52        0 (1260/ 326-1298/ 386)  80.5 E03-W00 S45-S39
  

    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