Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1511 (N42E07 L=351)  1571    313        4 ( 682/1450- 911/1608) -46.8 E13-W02 N35-N51
      2 1201 1352 (N22W18 L= 16)   289    116        0 (1184/1305-1225/1403)  71.6 W17-W20 N19-N25
      3  245 1294 (N36E59 L=298)  1366    155        3 ( 174/1242- 309/1332)  -9.5 E74-E50 N32-N38
      4  439 1157 (N25E37 L=321)   951    172        3 ( 365/1092- 482/1211) -45.0 E43-E33 N20-N29
      5  193  940 (N15E64 L=294)   760    121        2 ( 175/ 882- 213/1001) -90.0 E66-E61 N11-N18
      6  598  966 (N11E27 L=330)   256     84        1 ( 555/ 947- 635/ 984)  -9.5 E30-E24 N10-N11
      7  262  778 (N04E59 L=299)   405     72        3 ( 237/ 746- 290/ 806) -90.0 E61-E55 N02-N06
      8  827  684 (S09E19 L=339)   812    298        1 ( 717/ 602- 902/ 768)  25.0 E27-E12 S12-S05
      9  536  596 (S10E40 L=317)   292     31        0 ( 513/ 577- 559/ 617) -18.4 E41-E38 S11-S08
     10 1453  563 (S31W19 L= 16)   305     92        0 (1412/ 549-1497/ 581)   0.0 W15-W22 S32-S29
     11 1297  434 (S36W05 L=  3)   550     79        2 (1264/ 404-1333/ 467)  26.6 W02-W08 S38-S34
     12  952  316 (S37E23 L=334)   378     48        0 ( 933/ 286- 974/ 343)  38.7 E25-E20 S39-S35
  

    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