Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  617 1798 (N60E35 L=190)   676     80        2 ( 577/1772- 658/1814)   4.8 E41-E28 N57-N62
      2 1615 1603 (N29W54 L=279)  1565    212        1 (1529/1520-1672/1660) -65.6 W49-W58 N23-N34
      3  919 1605 (N37E01 L=224)   451     72        0 ( 889/1581- 961/1626) -36.9 E02-W02 N35-N38
      4 1295 1521 (N27W24 L=249)   265     75        0 (1260/1501-1326/1537) -45.0 W21-W26 N25-N28
      5 1780 1332 (N10W58 L=283)   634    122        1 (1756/1272-1801/1386)  81.9 W56-W60 N06-N13
      6  304  904 (S01E50 L=175)   848    157        0 ( 235/ 841- 370/ 947) -35.0 E56-E44 S05-N02
      7  273  638 (S16E61 L=163)  1245    117        5 ( 229/ 599- 317/ 682) -33.7 E66-E56 S18-S13
      8  381  631 (S18E50 L=174)   361     91        1 ( 350/ 598- 409/ 673)  26.6 E54-E46 S19-S15
      9 1287  528 (S35W13 L=238)   325    106        0 (1270/ 476-1312/ 580) -76.0 W11-W15 S38-S30
     10  466  532 (S25E47 L=178)   853    221        1 ( 411/ 452- 556/ 574)  17.5 E56-E37 S29-S23
     11 1190  298 (S52W04 L=229)   593     89        2 (1146/ 275-1238/ 320)  12.5 W00-W09 S53-S50
     12  955  190 (S59E25 L=200)  1149    102        3 ( 909/ 153- 996/ 232)  20.6 E32-E16 S63-S55
     13 1267  173 (S68W16 L=241)   474     76        0 (1234/ 160-1311/ 190)  -4.4 W10-W23 S70-S65
  

    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