Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1605 (N41W04 L= 88)   460     76        1 (1219/1585-1306/1621)  15.9 W00-W07 N38-N42
      2  511 1550 (N21E48 L= 35)  1115    283        5 ( 399/1506- 617/1597)  26.6 E57-E40 N16-N25
      3  979 1535 (N29E14 L= 69)   240    134        0 ( 931/1508-1033/1580) -12.5 E18-E09 N28-N31
      4 1628 1359 (N32W37 L=120)   402     64        1 (1605/1328-1651/1392) -53.1 W34-W38 N29-N33
      5 1360 1261 (N19W16 L=100)   289    106        2 (1324/1237-1400/1287) -18.4 W13-W19 N17-N20
      6  201  977 (S19E61 L= 22)  2299    355        2 ( 118/ 847- 274/1126)  36.5 E82-E53 S28-S08
      7 1798 1036 (N15W54 L=138)   399     74        1 (1772/1007-1814/1073)  90.0 W51-W56 N12-N17
      8 1102 1008 (S01W05 L= 88)   576    180        1 (1010/ 987-1174/1030)  18.4 W00-W09 S03-N01
      9  818  821 (S17E09 L= 75)   531    256        0 ( 697/ 790- 932/ 846)  25.0 E16-E01 S21-S13
     10 1429  682 (S14W33 L=116)   749    272        2 (1318/ 596-1527/ 753)  52.4 W27-W37 S20-S07
     11  443  675 (S34E35 L= 48)   799    219        4 ( 371/ 603- 532/ 736)  38.7 E44-E27 S40-S28
     12 1766  646 (S09W63 L=146)  1384    171        3 (1727/ 567-1797/ 722) -79.7 W60-W66 S14-S03
     13 1699  593 (S13W59 L=142)   430     92        1 (1676/ 544-1723/ 629) -80.5 W57-W60 S16-S10
     14  346  564 (S45E50 L= 34)   413     35        0 ( 321/ 550- 372/ 578)   0.0 E54-E45 S45-S42
     15 1631  535 (S18W55 L=138)   352     96        1 (1589/ 500-1664/ 572) -15.3 W49-W60 S19-S16
     16 1518  489 (S24W48 L=131)   315     73        0 (1484/ 470-1544/ 519)  53.1 W45-W49 S25-S21
  

    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