Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1582 1524 (N42W33 L=302)   328    100        0 (1557/1495-1617/1538)  39.8 W30-W36 N39-N44
      2 1311 1244 (N17W13 L=283)  3443    738        5 (1031/ 952-1537/1525) -31.4 E09-W33 N03-N29
      3  297 1343 (N05E58 L=211)   531    101        1 ( 259/1322- 328/1375)  -1.3 E61-E54 N03-N06
      4  201 1050 (S14E60 L=209)   576     50        1 ( 177/1020- 222/1075)  76.0 E62-E57 S16-S12
      5  128 1005 (S18E69 L=200)  1252    129        1 ( 109/ 964- 154/1054)  80.5 E71-E65 S20-S14
      6 1568  863 (S00W37 L=306)   244     54        0 (1542/ 840-1594/ 884)   0.0 W35-W39 S01-N00
      7  537  779 (S25E27 L=242)   385    113        1 ( 487/ 762- 592/ 795)   7.1 E31-E23 S25-S23
      8 1311  766 (S11W22 L=291)   740    202        1 (1209/ 724-1408/ 795)  36.0 W16-W28 S15-S07
      9 1045  702 (S20W08 L=277)   773    177        0 ( 972/ 679-1144/ 743)   5.2 W02-W13 S21-S18
     10 1090  429 (S37W19 L=288)   248     44        1 (1073/ 408-1106/ 455) -45.0 W17-W20 S38-S35
     11 1462  312 (S37W53 L=322)   985    158        1 (1433/ 283-1493/ 364) -29.1 W47-W56 S39-S33
     12 1274  311 (S41W38 L=307)  2160    400        2 (1095/ 284-1424/ 345)  12.1 W23-W51 S44-S36
  

    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