Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1263 1733 (N54W14 L= 83)   350     65        1 (1227/1719-1298/1747)  14.0 W09-W17 N52-N55
      2  558 1490 (N26E39 L= 30)  1481    344        5 ( 468/1403- 607/1619) -78.7 E45-E35 N19-N34
      3 1486 1373 (N28W29 L= 98)  2339    586        2 (1286/1179-1649/1514) -35.1 W14-W41 N16-N35
      4 1292 1211 (N15W15 L= 84)   397    146        0 (1262/1153-1337/1279) -54.5 W12-W18 N12-N19
      5 1484 1188 (N16W29 L= 98)   288    121        1 (1465/1136-1495/1245) -81.9 W27-W29 N12-N19
      6 1691 1204 (N18W46 L=115)   391    108        0 (1669/1172-1729/1244) -36.9 W43-W49 N16-N20
      7 1026  931 (S04W01 L= 70)   220     91        0 ( 987/ 908-1053/ 970)  63.4 E01-W02 S05-S01
      8  392  762 (S21E42 L= 27)   345     94        1 ( 373/ 725- 422/ 810) -90.0 E43-E39 S23-S18
      9 1394  650 (S18W29 L= 98)   768    206        1 (1313/ 590-1472/ 706)  45.0 W23-W32 S22-S14
     10 1670  593 (S20W53 L=123)  1092    186        1 (1633/ 516-1716/ 674) -84.8 W50-W56 S25-S14
     11 1725  593 (S19W60 L=129)   551     72        1 (1698/ 555-1753/ 628)  80.5 W57-W62 S22-S16
     12  483  455 (S43E41 L= 28)   437     68        0 ( 454/ 428- 511/ 485)  33.7 E46-E37 S45-S39
  

    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