Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1545 (N27E36 L=226)   395     89        0 ( 533/1500- 574/1582) -78.7 E37-E34 N24-N29
      2 1690 1429 (N25W46 L=307)  1767    312        2 (1578/1298-1773/1535) -47.1 W37-W51 N16-N31
      3  444 1475 (N22E42 L=219)   381     95        0 ( 407/1449- 475/1496) -26.6 E45-E39 N20-N23
      4 1369 1421 (N22W21 L=282)   230     50        0 (1340/1405-1397/1436)  14.0 W18-W22 N20-N22
      5  257 1337 (N14E56 L=205)   382     51        1 ( 245/1310- 273/1367) -63.4 E58-E54 N11-N15
      6 1623 1257 (N13W38 L=299)  1329    220        1 (1548/1188-1718/1348) -36.9 W32-W45 N09-N18
      7  733 1298 (N11E19 L=242)   230     62        1 ( 705/1267- 753/1330) -63.4 E21-E17 N09-N13
      8 1206 1264 (N11W10 L=271)   294    123        0 (1145/1240-1272/1287)  14.0 W06-W14 N09-N12
      9 1817 1156 (N09W55 L=316)   333     67        0 (1805/1123-1830/1187)  76.0 W53-W56 N06-N10
     10 1859 1107 (N06W60 L=321)  1157    222        4 (1816/1064-1903/1154) -45.0 W54-W65 N03-N09
     11  154 1054 (S03E65 L=196)  1144    215        2 ( 129/ 952- 176/1141)  73.3 E68-E61 S09-N01
     12  376 1011 (S06E42 L=219)   945    198        1 ( 342/ 920- 442/1074)  60.9 E45-E37 S11-S02
     13  535  960 (S09E31 L=231)   334    104        2 ( 500/ 923- 582/ 995)  33.7 E33-E27 S11-S07
     14  728  799 (S19E18 L=243)  1064    293        0 ( 630/ 693- 866/ 872)  34.5 E25-E09 S25-S13
     15 1714  789 (S14W48 L=310)   485     98        1 (1673/ 758-1744/ 833) -53.1 W44-W51 S15-S11
     16 1537  609 (S27W38 L=299)  1341    253        1 (1428/ 557-1613/ 716)  42.3 W30-W41 S30-S19
     17  596  545 (S35E30 L=231)   692    221        1 ( 540/ 474- 660/ 604)  39.3 E35-E23 S40-S31
     18 1766  566 (S27W60 L=321)   828    110        0 (1716/ 529-1803/ 591)  39.8 W56-W63 S30-S25
  

    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