Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1332 1675 (N43E02 L=198)   661    189        0 (1264/1596-1386/1744) -28.3 E09-W04 N39-N46
      2  489 1364 (S01E42 L=157)  1228    264        3 ( 444/1245- 554/1466) -37.6 E49-E34 S07-N04
      3 1405 1263 (N18W15 L=214)   727    170        0 (1335/1197-1458/1310)  -5.7 W09-W19 N15-N19
      4  979  811 (S19W04 L=203)  1172    454        0 ( 780/ 731-1179/ 852)  45.0 E05-W13 S29-S11
      5 1572  829 (S01W38 L=237)   391     78        0 (1532/ 807-1617/ 850)   8.1 W34-W41 S02-N00
      6  686  814 (S27E14 L=185)   205     70        0 ( 658/ 786- 712/ 845)  68.2 E15-E13 S29-S24
      7  541  690 (S39E22 L=177)  1172    292        3 ( 497/ 566- 611/ 789) -82.4 E26-E17 S46-S31
      8 1077  649 (S25W15 L=214)   616    158        0 (1011/ 612-1161/ 678)  53.1 W11-W18 S29-S21
      9 1415  609 (S17W36 L=236)   370     72        1 (1379/ 584-1450/ 631)  53.1 W34-W38 S18-S14
     10 1199  540 (S28W27 L=226)   771    199        0 (1105/ 517-1298/ 566)  35.0 W21-W31 S31-S23
     11  458  419 (S63E33 L=167)  1153    157        1 ( 395/ 387- 527/ 469)  -5.4 E42-E21 S65-S59
     12 1421  407 (S27W49 L=248)   802    185        0 (1348/ 377-1498/ 446)   7.1 W41-W57 S28-S25
  

    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