Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1075 1432 (N30W05 L=300)   684    229        1 (1040/1333-1134/1495)  70.0 W02-W10 N23-N34
      2 1657 1421 (N26W51 L=346)   649    128        1 (1614/1368-1688/1476) -74.1 W49-W53 N21-N29
      3  270 1349 (N25E61 L=234)   758     98        2 ( 240/1306- 299/1389)  51.3 E63-E58 N22-N28
      4 1154 1319 (N22W10 L=304)   329     91        1 (1107/1306-1204/1339)  -8.1 W06-W13 N21-N23
      5  334 1102 (N10E48 L=246)   729    185        2 ( 304/1034- 363/1175)  83.7 E51-E45 N05-N14
      6 1749  856 (S10W51 L=346)   495     34        1 (1722/ 834-1772/ 879) -21.8 W48-W53 S11-S08
      7 1003  820 (S08E02 L=293)   375     70        1 ( 980/ 779-1022/ 858) -78.7 E03-W00 S10-S05
      8  291  671 (S17E58 L=237)   395     64        1 ( 264/ 646- 318/ 700) -21.8 E60-E55 S18-S15
      9  441  647 (S18E43 L=251)   319     37        0 ( 422/ 624- 457/ 670)  36.9 E45-E41 S20-S16
     10 1657  606 (S26W47 L=341)   747    119        0 (1634/ 576-1696/ 654) -36.9 W43-W51 S28-S22
     11 1663  492 (S35W53 L=348)   594     59        1 (1628/ 471-1693/ 518)   7.1 W48-W56 S36-S33
  

    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