Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1214 1205 (N13W08 L=103)   507    145        0 (1167/1159-1278/1266) -24.0 W03-W12 N11-N15
      2  591  983 (S12E26 L= 70)   418    119        0 ( 576/ 925- 609/1054) -66.8 E27-E24 S15-S08
      3 1177  967 (S02W10 L=106)   196    104        0 (1151/ 912-1204/1014) -38.7 W07-W12 S04-N00
      4 1791  918 (N08W55 L=150)   631    151        1 (1747/ 869-1825/ 951) -15.3 W49-W60 N05-N09
      5 1368  885 (S03W23 L=119)   302     68        0 (1348/ 844-1390/ 922)  90.0 W22-W24 S05-N00
      6  889  660 (S26E01 L= 95)  1368    538        1 ( 650/ 548-1132/ 748)  38.2 E16-W12 S38-S16
      7 1376  617 (S18W31 L=127)   284    101        0 (1335/ 568-1420/ 641)  63.4 W29-W33 S22-S16
      8 1539  561 (S18W45 L=141)   402     89        1 (1511/ 530-1562/ 600)  90.0 W43-W46 S20-S15
      9  904  487 (S37W04 L=100)   274     67        0 ( 889/ 462- 925/ 519) -53.1 W03-W06 S38-S34
     10 1555  461 (S24W51 L=147)   388    144        0 (1499/ 449-1605/ 471)  11.3 W46-W56 S25-S22
     11 1381  268 (S41W51 L=146)  1375    251        2 (1284/ 244-1478/ 295)  28.3 W44-W58 S45-S38
  

    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