Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  460 1665 (N17E58 L=278)   456    119        0 ( 444/1633- 477/1702) -51.3 E60-E56 N15-N20
      2  497 1612 (N14E51 L=284)   457    112        0 ( 472/1584- 530/1645) -31.0 E53-E47 N12-N16
      3  490 1534 (N09E48 L=288)   323    118        0 ( 470/1493- 512/1569) -31.0 E50-E44 N07-N11
      4 1103 1276 (N09E02 L=334)   449    140        2 (1051/1205-1140/1332) -35.5 E06-W01 N05-N10
      5  417 1218 (S10E42 L=294)   759    179        1 ( 346/1184- 478/1248)  21.8 E47-E37 S12-S08
      6 1497 1123 (N11W24 L=  0)   414    126        1 (1478/1052-1513/1191) -63.4 W22-W26 N07-N15
      7 1342 1026 (N02W18 L=353)   598    235        2 (1259/ 910-1387/1121) -30.3 W10-W22 S03-N04
      8  611  839 (S28E20 L=316)   628    254        1 ( 506/ 814- 742/ 861)  13.2 E28-E11 S30-S25
      9  412  785 (S36E35 L=301)   300     87        1 ( 374/ 766- 449/ 809)  23.2 E38-E31 S37-S33
     10 1579  767 (S05W40 L= 16)   247     92        0 (1539/ 732-1601/ 785)   0.0 W36-W42 S06-S04
     11  845  567 (S38W03 L=339)   287     88        0 ( 802/ 537- 882/ 600)  66.8 W00-W04 S41-S34
     12 1522  402 (S26W54 L= 30)   451     88        1 (1496/ 386-1559/ 428)   8.1 W50-W57 S26-S24
  

    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