Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1269 1708 (N47W34 L=134)  1472    305        4 (1240/1599-1321/1783)  39.1 W28-W44 N38-N52
      2  937 1501 (N37W01 L=100)   295     46        1 ( 904/1477- 965/1517) -33.7 E01-W02 N35-N38
      3 1001 1425 (N31W04 L=103)   222     58        0 ( 983/1399-1025/1454) -63.4 W02-W05 N29-N33
      4  735 1371 (N31E16 L= 84)   406    140        0 ( 700/1295- 764/1425) -66.0 E18-E13 N25-N34
      5 1093 1286 (N21W08 L=107)   236     74        0 (1057/1255-1129/1320) -59.0 W05-W09 N18-N23
      6 1577 1214 (N09W39 L=138)   283     79        1 (1540/1181-1599/1247) -59.0 W36-W40 N07-N12
      7 1249 1177 (N12W16 L=115)   355    103        0 (1204/1137-1282/1230) -63.4 W13-W17 N09-N15
      8 1094 1161 (N13W06 L=105)   264     95        2 (1043/1128-1139/1187)  18.4 W02-W08 N11-N14
      9  240 1036 (N13E57 L= 43)   590    112        1 ( 210/1007- 277/1077) -33.7 E59-E52 N11-N15
     10  943  699 (S13E09 L= 90)   219     70        0 ( 929/ 657- 955/ 735)  68.2 E10-E07 S15-S10
     11 1452  560 (S30W24 L=123)  1492    275        1 (1325/ 481-1573/ 667)  25.3 W13-W33 S33-S24
     12  637  477 (S25E35 L= 65)   869    186        0 ( 564/ 436- 702/ 529)  21.0 E41-E28 S27-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