Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1132 1872 (N53E29 L=219)   400     36        0 (1106/1858-1158/1888)  76.0 E30-E27 N50-N54
      2 1289 1794 (N52E10 L=238)   746     81        1 (1248/1763-1335/1826)  -6.3 E14-E05 N49-N54
      3  397 1625 (N13E63 L=185)   734     54        0 ( 368/1599- 428/1649)  56.3 E66-E60 N10-N14
      4 1561 1620 (N48W20 L=268)   388     51        0 (1535/1595-1589/1643) -21.8 W17-W22 N47-N49
      5 1661 1442 (N38W31 L=279)   435     36        0 (1638/1421-1684/1467) -11.3 W28-W33 N36-N40
      6  267 1436 (S01E65 L=183)   571     87        0 ( 236/1394- 286/1462) -15.9 E68-E60 S03-N00
      7  218 1337 (S07E66 L=182)  1293    102        2 ( 183/1288- 245/1376)  90.0 E72-E62 S10-S04
      8  200 1199 (S16E63 L=185)  2118    317        4 (  88/1131- 352/1257)  13.7 E84-E47 S20-S10
      9  643 1185 (S07E26 L=222)   627    248        0 ( 602/1098- 734/1242)  78.7 E27-E20 S12-S02
     10 1313 1194 (N11W11 L=259)   437     94        0 (1268/1148-1351/1236) -31.0 W07-W14 N08-N12
     11  967 1090 (S04E05 L=243)   276     42        0 ( 945/1064- 988/1120)  90.0 E05-E03 S05-S02
     12  436  954 (S25E35 L=213)   677    164        0 ( 401/ 878- 491/1017)  68.2 E37-E31 S30-S20
     13 1165  798 (S14W14 L=262)   251     94        1 (1132/ 760-1198/ 829) -36.9 W11-W15 S16-S12
     14 1376  669 (S15W31 L=279)   297     72        0 (1340/ 646-1398/ 701)  78.7 W29-W32 S17-S12
     15  525  544 (S49E22 L=226)  1283    250        2 ( 403/ 513- 637/ 601)   6.8 E35-E10 S51-S45
     16 1195  528 (S28W26 L=274)   515    118        0 (1147/ 507-1260/ 559)  14.0 W22-W30 S29-S26
     17 1530  421 (S23W55 L=303)  1439    238        4 (1427/ 389-1637/ 454)  26.6 W46-W63 S27-S18
     18 1125  253 (S47W42 L=290)  1565    285        3 (1010/ 186-1258/ 321)   7.9 W28-W59 S49-S44
  

    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