Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1378 1631 (N44W30 L=130)  1110    222        3 (1327/1566-1419/1690)  33.7 W24-W36 N38-N48
      2  904 1469 (N32E09 L= 91)   341    101        0 ( 889/1419- 920/1519) -81.9 E09-E07 N28-N35
      3 1713 1458 (N30W55 L=154)   441     98        2 (1691/1434-1742/1480) -26.6 W52-W57 N28-N31
      4  544 1311 (N21E33 L= 66)   433     97        1 ( 528/1261- 565/1362)  90.0 E34-E31 N17-N23
      5 1529 1314 (N21W35 L=134)   449    116        1 (1478/1279-1584/1351) -35.5 W31-W38 N18-N23
      6 1530 1224 (N15W34 L=133)   495    164        1 (1462/1205-1595/1254)   0.0 W28-W39 N14-N17
      7 1072 1160 (N12W03 L=102)   318    105        0 (1052/1113-1103/1216) -66.8 W01-W04 N08-N15
      8 1516 1111 (N08W32 L=132)   481    235        0 (1426/1076-1612/1188) -28.3 W26-W39 N06-N13
      9 1251 1111 (N09W14 L=114)   291    155        2 (1177/1059-1321/1155) -33.7 W09-W18 N05-N11
     10  468  997 (N01E37 L= 63)   373    138        1 ( 416/ 966- 523/1036)  26.6 E40-E32 S01-N03
     11 1347  777 (S12W21 L=120)   228     80        0 (1321/ 757-1386/ 802) -36.9 W19-W23 S13-S10
     12  760  745 (S14E17 L= 83)   247     68        1 ( 743/ 718- 782/ 776)  56.3 E17-E15 S15-S12
     13  372  710 (S18E46 L= 53)   416    110        1 ( 352/ 663- 389/ 761)  63.4 E48-E44 S20-S14
     14 1559  628 (S22W38 L=137)   587    107        3 (1519/ 594-1599/ 663)  18.4 W34-W40 S24-S19
     15  688  628 (S22E23 L= 77)   627    126        1 ( 629/ 604- 751/ 659) -14.0 E26-E18 S23-S20
     16  495  608 (S24E38 L= 62)   778    225        4 ( 416/ 578- 574/ 657)  23.2 E44-E30 S26-S20
     17  698  377 (S41E27 L= 73)   403     89        1 ( 669/ 355- 724/ 411)  38.7 E30-E24 S43-S38
     18  963  353 (S43E05 L= 95)   922    141        3 ( 898/ 322-1004/ 408)  45.0 E10-E01 S45-S37
  

    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