Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  684 1504 (N31E28 L=354)   467    188        2 ( 617/1441- 749/1531) -22.6 E33-E21 N27-N33
      2 1572 1495 (N35W42 L= 64)   323     18        0 (1553/1471-1587/1517) -71.6 W39-W43 N33-N36
      3 1433 1475 (N33W29 L= 51)   414     80        0 (1393/1446-1467/1509) -31.0 W25-W31 N30-N35
      4 1054 1437 (N29W00 L= 22)   487    233        0 (1003/1373-1100/1478) -39.8 E03-W03 N24-N31
      5  888 1225 (N14E09 L= 13)   242     85        0 ( 851/1208- 928/1246) -21.8 E11-E06 N13-N15
      6  330 1163 (N07E49 L=333)   331     41        0 ( 310/1137- 349/1188)  56.3 E50-E47 N05-N08
      7  383  925 (S07E43 L=339)   789    234        2 ( 356/ 839- 417/1023) -90.0 E45-E40 S12-N00
      8  266  818 (S15E55 L=327)  2578    419        7 ( 179/ 687- 309/ 990)  48.0 E72-E49 S23-S03
      9  644  693 (S20E24 L=358)  1050    255        5 ( 577/ 598- 709/ 797) -69.0 E28-E19 S26-S13
     10  430  754 (S17E40 L=342)   272     62        2 ( 401/ 735- 467/ 781)  18.4 E42-E36 S18-S15
     11  710  737 (S17E19 L=  3)   358     54        2 ( 688/ 703- 740/ 766)  78.7 E20-E17 S19-S14
     12 1837  728 (S14W66 L= 88)   673     75        1 (1795/ 711-1872/ 744)  10.3 W60-W72 S14-S12
     13  260  605 (S29E64 L=318)  1529    170        2 ( 233/ 540- 304/ 682)  84.3 E69-E59 S33-S23
     14  888  538 (S29E07 L= 15)  2834    623        8 ( 630/ 433-1128/ 669)  28.5 E28-W08 S38-S19
     15  544  616 (S26E32 L=350)   434    129        0 ( 504/ 577- 604/ 644)  29.1 E36-E27 S28-S23
     16  351  565 (S31E53 L=329)   539     87        1 ( 307/ 542- 391/ 583)   0.0 E59-E48 S33-S29
     17 1297  497 (S30W22 L= 45)   436    165        0 (1229/ 461-1367/ 530)  33.7 W17-W27 S33-S27
  

    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