Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1586 1603 (N28W59 L=246)   546    163        0 (1518/1561-1627/1635) -49.4 W53-W63 N24-N31
      2 1641 1507 (N21W57 L=244)   488     67        3 (1614/1481-1667/1541)  15.9 W53-W60 N20-N23
      3 1074 1355 (N26W10 L=197)  1033    263        1 ( 987/1256-1154/1451) -68.2 W05-W12 N18-N33
      4 1289 1382 (N23W24 L=212)   370    170        0 (1210/1364-1360/1403) -11.3 W19-W29 N21-N24
      5  651 1111 (N18E22 L=165)   249    115        0 ( 623/1068- 673/1168) -74.1 E23-E20 N15-N22
      6  870 1053 (N11E09 L=178)   243    149        2 ( 840/1011- 933/1089)  15.9 E11-E04 N08-N13
      7  997  827 (S04E05 L=182)   405     97        1 ( 955/ 792-1042/ 862) -56.3 E07-E02 S07-S01
      8 1126  672 (S16E00 L=187)   324    133        0 (1096/ 625-1145/ 744)  60.3 E03-W02 S18-S11
      9 1304  590 (S25W10 L=197)  1073    400        0 (1187/ 478-1455/ 716)  25.5 W00-W21 S30-S20
     10 1474  587 (S30W22 L=209)   932    322        1 (1371/ 468-1568/ 681)  32.0 W12-W29 S36-S25
     11  350  546 (S13E59 L=128)  1082    163        3 ( 321/ 483- 383/ 615)  83.7 E63-E56 S17-S08
     12  858  485 (S24E22 L=166)   285     61        0 ( 827/ 458- 884/ 513)  21.8 E24-E19 S25-S21
  

    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