Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  474 1668 (N34E48 L=132)  1076    260        0 ( 428/1626- 526/1720) -28.1 E56-E41 N30-N38
      2 1619 1416 (N25W38 L=219)   669    122        0 (1558/1388-1665/1435)   0.0 W33-W41 N22-N26
      3 1827 1287 (N19W56 L=237)   522     53        0 (1804/1254-1849/1322) -45.0 W53-W58 N17-N21
      4  395 1266 (N07E43 L=138)   684    119        0 ( 373/1217- 427/1318) -63.4 E44-E40 N03-N09
      5 1922 1116 (N10W67 L=248)  1399    180        2 (1888/1032-1956/1180) -41.6 W62-W74 N05-N13
      6 1303  996 (S04W17 L=197)  1022    313        1 (1246/ 859-1339/1149) -79.4 W12-W18 S12-N04
      7 1671  951 (S04W43 L=223)   378    106        0 (1622/ 928-1717/ 969)  14.0 W38-W46 S05-S02
      8  798  883 (S15E13 L=168)   204     93        0 ( 755/ 857- 854/ 897)  18.4 E15-E09 S16-S13
      9 1690  807 (S12W46 L=227)   340     96        1 (1651/ 777-1710/ 837) -11.3 W43-W48 S13-S09
     10  442  792 (S22E38 L=143)   349    101        1 ( 403/ 764- 489/ 816)  23.2 E40-E33 S23-S20
     11 1808  639 (S20W63 L=243)  2914    348        4 (1705/ 556-1898/ 707)  38.7 W54-W73 S26-S14
     12 1460  418 (S39W40 L=221)   793    208        1 (1375/ 398-1554/ 447)  26.6 W33-W46 S41-S35
     13 1319  282 (S51W37 L=217)   331     56        1 (1288/ 264-1341/ 304)  33.7 W33-W38 S52-S48
  

    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