Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1265 1767 (N56W15 L= 18)   404     87        2 (1222/1756-1317/1778)   9.5 W09-W21 N54-N58
      2  450 1665 (N37E61 L=302)  2126    156        4 ( 423/1593- 483/1739) -17.2 E81-E51 N31-N41
      3  633 1624 (N36E38 L=325)  2329    275        6 ( 539/1563- 760/1697)   0.0 E48-E25 N31-N41
      4  539 1447 (N23E38 L=325)   250     62        0 ( 527/1417- 554/1482) -76.0 E39-E37 N21-N25
      5  606 1401 (N21E32 L=331)   618    182        1 ( 535/1357- 652/1478)  66.0 E36-E29 N17-N26
      6 1207 1385 (N25W09 L= 12)   273     98        0 (1188/1341-1240/1420)  71.6 W07-W10 N22-N28
      7 1708 1299 (N24W48 L= 51)   301     77        0 (1692/1259-1724/1337) -78.7 W46-W49 N21-N26
      8 1006 1322 (N20E04 L=359)   244     93        1 ( 955/1307-1059/1331)  15.9 E07-W00 N18-N20
      9 1434 1195 (N15W25 L= 28)   241    102        0 (1389/1185-1484/1215)  -8.1 W21-W28 N14-N16
     10 1764 1069 (N10W51 L= 55)   634    157        1 (1736/ 988-1788/1133) -58.0 W48-W54 N05-N13
     11 1418  991 (N03W25 L= 28)   497    187        0 (1365/ 930-1491/1021)  23.2 W21-W29 S01-N04
     12  237  910 (S13E56 L=307)   878    176        6 ( 206/ 834- 259/ 971) -63.4 E59-E53 S17-S08
     13 1752  722 (S11W56 L= 59)  3456    501        3 (1620/ 589-1870/ 944)  54.0 W47-W64 S20-N02
     14 1522  875 (S03W33 L= 36)   259     39        1 (1502/ 856-1551/ 897)  14.0 W31-W35 S04-S02
     15  159  862 (S17E68 L=295)   531     56        0 ( 134/ 839- 195/ 879)  -5.7 E72-E62 S17-S15
     16  427  594 (S32E42 L=322)  1165    230        2 ( 356/ 512- 514/ 657)  29.7 E53-E32 S38-S26
  

    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