Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1484 1591 (N23W49 L=179)   473    138        1 (1422/1551-1541/1618) -60.3 W45-W51 N19-N26
      2 1331 1459 (N21W31 L=162)  1184    487        0 (1123/1407-1530/1525) -36.0 W21-W44 N14-N30
      3 1790 1322 (S01W60 L=190)   554     74        1 (1752/1296-1821/1349) -36.9 W57-W63 S02-N01
      4 1093 1229 (N15W09 L=140)   352    148        0 (1029/1216-1167/1239) -20.6 W05-W13 N12-N16
      5  779 1177 (N20E10 L=120)   378    113        1 ( 722/1156- 840/1193) -23.2 E13-E06 N18-N22
      6 1891 1029 (S20W62 L=192)  1139    170        5 (1866/ 963-1913/1119)  66.0 W57-W65 S23-S14
      7  291 1086 (N27E49 L= 82)   401     83        2 ( 252/1063- 330/1107) -20.6 E52-E44 N25-N28
      8  181 1063 (N27E62 L= 68)   567     96        2 ( 157/1025- 217/1103) -39.8 E66-E57 N24-N29
      9  500  994 (N17E32 L= 98)   232    102        0 ( 453/ 983- 553/1003) -14.0 E36-E28 N15-N17
     10 1273  624 (S24W04 L=134)   910    233        2 (1241/ 526-1290/ 751)  61.7 W00-W07 S29-S16
     11 1227  464 (S33E04 L=126)   272     46        0 (1203/ 448-1259/ 480) -14.0 E05-E01 S34-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