Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  681 1522 (N18E30 L=246)   463     97        0 ( 638/1499- 730/1543)  21.8 E32-E27 N16-N19
      2 1710 1442 (N33W42 L=318)  1278    202        2 (1641/1346-1771/1524) -45.0 W35-W48 N27-N37
      3  365 1370 (N05E49 L=227)   427     89        1 ( 345/1323- 377/1411) -59.0 E51-E47 N02-N07
      4  150 1161 (S09E67 L=209)  1203    161        2 ( 122/1097- 171/1247)  90.0 E70-E64 S12-S03
      5 1846 1157 (N18W55 L=331)   833    154        2 (1772/1126-1916/1179)  12.5 W46-W64 N15-N21
      6 1061 1126 (S00W00 L=276)   208     49        1 (1026/1109-1088/1142)   0.0 E01-W01 N00-N01
      7 1894  989 (N09W62 L=338)  4271    244       13 (1833/ 875-1957/1097) -57.5 W54-W72 N01-N14
      8 1013  712 (S25W05 L=281)   454    143        0 ( 965/ 651-1068/ 771)  66.0 W02-W07 S29-S20
      9 1177  631 (S27W17 L=293)  1500    347        2 (1040/ 564-1349/ 681)  29.1 W09-W28 S32-S22
     10 1672  447 (S26W64 L=340)  3780    320        6 (1605/ 371-1767/ 562) -36.3 W54-W82 S30-S19
     11 1145  498 (S36W19 L=295)  1426    336        1 (1001/ 467-1306/ 545)  10.8 W08-W29 S38-S31
     12 1112  365 (S45W23 L=299)  3124    722        5 ( 829/ 260-1446/ 470)   1.8 E02-W60 S47-S41
  

    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