Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1686 1595 (N22W69 L= 17)   692     76        0 (1665/1557-1708/1633)  26.6 W66-W76 N19-N23
      2  750 1406 (N36E10 L=298)   898    255        0 ( 625/1379- 858/1431) -16.4 E19-E02 N32-N38
      3 1088 1362 (N25W12 L=320)   559    200        0 (1001/1307-1160/1391) -41.2 W06-W15 N20-N28
      4 1048 1192 (N16W05 L=313)   333    193        0 (1016/1136-1104/1230)  59.0 W02-W07 N13-N18
      5  779 1060 (N14E14 L=294)   205    109        0 ( 736/1029- 824/1089) -45.0 E16-E11 N11-N16
      6 1626  987 (S09W37 L=344)   538    182        0 (1566/ 920-1660/1056)  36.9 W30-W39 S11-S05
      7  434 1013 (N17E38 L=269)   263     62        0 ( 415/ 992- 465/1038) -33.7 E39-E35 N15-N18
      8  212  899 (N13E62 L=246)  1815    167        4 ( 180/ 836- 251/ 963)  60.3 E66-E56 N09-N16
      9  433  873 (N09E41 L=267)   465     91        1 ( 414/ 829- 458/ 918) -78.7 E42-E39 N06-N11
     10 1268  886 (S06W11 L=319)   274    181        0 (1201/ 844-1324/ 914) -14.0 W06-W14 S09-S04
     11  182  705 (N01E74 L=234)   745     66        2 ( 155/ 687- 214/ 726) -14.0 E81-E68 N00-N02
  

    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