Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1173 1471 (N24W10 L=357)   472    207        0 (1118/1406-1218/1534)  63.4 W06-W13 N20-N28
      2  472 1370 (N19E37 L=310)   724    158        1 ( 437/1294- 503/1442)  90.0 E40-E35 N14-N23
      3 1743 1294 (N14W51 L= 38)   396     57        0 (1722/1262-1761/1325) -56.3 W49-W52 N11-N15
      4  977 1171 (N06E03 L=344)   773    219        2 ( 916/1088-1029/1274)  79.7 E06-W00 N00-N11
      5  304  971 (S05E49 L=298)   274     48        0 ( 284/ 950- 328/ 996) -33.7 E50-E46 S05-S03
      6  677  882 (S11E22 L=325)   250     73        0 ( 649/ 853- 700/ 910) -45.0 E23-E20 S12-S09
      7  767  688 (S23E17 L=330)   689    258        1 ( 717/ 607- 794/ 782)  61.4 E21-E15 S28-S17
      8  202  715 (S19E66 L=281)   867     70        2 ( 173/ 682- 226/ 757) -90.0 E71-E62 S21-S16
      9  570  487 (S36E37 L=310)   901    145        1 ( 497/ 449- 638/ 517)  11.3 E45-E30 S38-S34
     10  795  431 (S41E19 L=328)   247     55        0 ( 765/ 418- 832/ 447)  18.4 E21-E15 S41-S39
     11 1485  340 (S48W45 L= 32)  1689    159        3 (1423/ 301-1560/ 400) -15.6 W35-W60 S51-S43
     12 1346  339 (S49W30 L= 17)   562     89        1 (1311/ 313-1377/ 377)  68.2 W25-W32 S50-S45
     13  783  229 (S58E30 L=317)  1214    186        1 ( 708/ 188- 863/ 263)   9.1 E43-E18 S62-S55
  

    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