Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  953 1702 (N51W19 L=291)   482     73        1 ( 928/1656- 976/1740)  36.9 W15-W23 N47-N53
      2 1368 1659 (N34W46 L=318)   315     38        0 (1341/1643-1400/1671) -26.6 W43-W48 N33-N35
      3 1644 1569 (N20W62 L=334)  1945    158        6 (1578/1499-1696/1611)  76.0 W57-W70 N16-N24
      4  416 1325 (N39E38 L=234)  1154    179        0 ( 333/1278- 487/1356) -26.6 E47-E30 N34-N42
      5  897 1257 (N24E02 L=270)  1193    401        2 ( 724/1224-1078/1284) -23.2 E13-W08 N17-N27
      6 1655 1108 (S05W41 L=313)  1061    192        0 (1621/1018-1683/1188)  41.6 W36-W45 S08-N00
      7 1209 1131 (N09W13 L=285)   374    132        0 (1159/1102-1252/1177)  15.9 W09-W16 N07-N10
      8  714 1126 (N20E16 L=255)   239     89        1 ( 672/1100- 747/1155)   0.0 E19-E14 N18-N21
      9  400 1052 (N22E41 L=231)  1045    186        2 ( 336/ 987- 469/1109) -42.0 E46-E35 N16-N26
     10 1500  532 (S36W19 L=291)  6752   1124       16 (1114/ 207-1759/ 914)  26.9 E20-W45 S52-S18
     11 1092  570 (S21E07 L=264)   515    193        0 (1035/ 509-1120/ 646)  41.2 E12-E04 S23-S16
     12  504  396 (S20E55 L=217)   418     56        0 ( 487/ 362- 517/ 424)  59.0 E57-E52 S22-S17
     13  945  330 (S34E26 L=246)   732    157        0 ( 890/ 291-1006/ 383)  18.4 E31-E19 S36-S31
  

    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